{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "194c31ee-ba2b-4cc7-995b-627255ccd35c",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-03-20T10:31:21.446027Z",
     "iopub.status.busy": "2022-03-20T10:31:21.445599Z",
     "iopub.status.idle": "2022-03-20T10:31:21.948808Z",
     "shell.execute_reply": "2022-03-20T10:31:21.948053Z",
     "shell.execute_reply.started": "2022-03-20T10:31:21.445916Z"
    },
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7f0991d54d10>]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVYAAAE/CAYAAAAdVOsRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XmcVNWd9/HPl2ZHdhAREFBAgktEiZKYREdcwCRiYhIxMWLCyDNuo0km45I86ugk0ehoNCGZGGVcYiS4TGxHjEGFOMkjCAREAYEGBRpRoFm12br5PX+c01DdVNHV9K2u6u7f+/XqV90699x7z+26/etb92wyM5xzziWnRb4L4JxzTY0HVuecS5gHVuecS5gHVuecS5gHVuecS5gHVuecS5gH1jyRtEjSmfkuh3MuefJ2rM45lyy/Y3XOuYR5YM0TSe9JOlvSqZJel7RF0jpJv5TUOiWfSfpnSSslbZR0t6QWcd0xkl6VVBbXPSGpS41j/IukhZK2SvqDpLb5OF/nmhMPrPlXCXwX6AF8GhgFXFUjz5eBEcDJwFjgOzFdwE+BI4FPAP2A22ps+3VgNDAQOBG4POHyO+dq8MCaZ2Y2z8xmmVmFmb0H/AY4o0a2u8xsk5mtBn4OXBK3LTGz6Wa2y8w2APem2fYBM3vfzDYBzwMn5fSEnHO0zHcBmjtJQwgBcQTQnvCZzKuRbU3K8irCHSqSegH3A58DOhL+UW6use0HKcvlVds653LH71jz79fAO8BgM+sE3Ez4ip+qX8ryUcD7cfkngAEnxG0vTbOtc66BeWDNv47ANuAjSUOBK9Pk+YGkrpL6AdcBf0jZ9iNgq6Q+wA8aosDOuYPzwJp//wJ8A9gO/Jb9QTPVc4THAwuAF4CHY/q/ESq0tsb0Z3NdWOdc7byDQJ5IWg1camav1ZLPCI8JShqmZM65+vI71jyQ1BPoCbyX56I453Igq8Aq6buxb/vbkp6U1FbSQEmzJZXEhuetY9428X1JXD8gZT83xfSlks5LSR8d00ok3Zj0SRYSSZ8ClgO/iM2nXD1ImixpvaS3M6yXpAfitbVQ0skNXUbX/NQaWGOlyD8DI8zseKAIGAfcBdxnZoMITXwmxE0mAJtj+n0xH5KGxe2OIzRY/5WkIklFwCRgDDAMuCTmbZLMbI6ZdTGzf80yv/wxwEE9QrieMhkDDI4/EwmtMJzLqWwfBbQE2klqSWhruQ44C3g6rn8UuDAuj43vietHSVJMnxIbs78LlACnxp8SM1tpZruBKTGvc7WKz6g3HSTLWOAxC2YBXST1bpjSueaq1sBqZmuBe4DVhIC6lVBDvcXMKmK2UqBPXO5DbNAe128Fuqem19gmU7pzSfDryzW4WnteSepK+K8/ENgCPMXBv3rljKSJhK9zdOjQ4ZShQ4fmoxguj+bNm7fRzHrmYt9+fTlI5hrLpkvr2cC7sS86kp4FTid8pWoZ70r7Amtj/rWEnkKl8dFBZ6AsJb1K6jaZ0qsxsweBBwFGjBhhc+fOzaL4rimRtKqOmxzsuqvGry8Hh3SNHSCbZ6yrgZGS2sdnpaOAxcAM4Ksxz3hCI3aA4vieuP5VC41li4FxsdXAQEJlwhvAHGBwbGXQmlDBVVzfE3MuKgYui60DRgJbzWxdvgvlmrZa71jNbLakp4G/AxXAfMJ/9ReAKZL+PaZV9QZ6GHhcUgmhUmFc3M8iSVMJQbkCuNrMKgEkXQO8RGhxMNnMFiV3iq4pk/QkcCbQQ1IpcCvQCsDM/hOYBpxPqCwtB76dn5K65qTR9rzyr2rNk6R5ZjYi18fx66v5SuIa855XzjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+szjmXMA+srtGTNFrSUkklkm5Ms/4oSTMkzZe0UNL5+Sinaz48sLpGTVIRMAkYAwwDLpE0rEa2HwFTzWw4YaqgXzVsKV1z44HVNXanAiVmttLMdgNTCNO1pzKgU1zuDLzfgOVzzVA20187V8j6AGtS3pcCp9XIcxvwZ0nXAh0IU7o7lzN+x+qag0uAR8ysL2HG1sclHXDtS5ooaa6kuRs2bGjwQrqmwwOra+zWAv1S3veNaakmAFMBzOx1oC3Qo+aOzOxBMxthZiN69uyZo+K65iCrwCqpi6SnJb0jaYmkT0vqJmm6pOXxtWvMK0kPxBrahZJOTtnP+Jh/uaTxKemnSHorbvOAJCV/qq6JmgMMljRQUmtC5VRxjTyrgVEAkj5BCKx+S+pyJts71vuBP5nZUOCTwBLgRuAVMxsMvBLfQ6idHRx/JgK/BpDUDbiV8PzrVODWqmAc81yRst3o+p2Way7MrAK4BniJcF1ONbNFkm6XdEHM9n3gCklvAk8Cl5uZ5afErjmotfJKUmfg88DlALHmdbekscCZMdujwEzgBkKN7GPxwp0V73Z7x7zTzWxT3O90YLSkmUAnM5sV0x8DLgReTOQMXZNnZtOAaTXSbklZXgyc3tDlcs1XNnesAwlfm/4rNrB+SFIHoJeZrYt5PgB6xeV0tbR9akkvTZPunHONUjaBtSVwMvDr2MD6Y/Z/7Qcg3p3m/KuV19o65xqDbAJrKVBqZrPj+6cJgfbD+BWf+Lo+rs9US3uw9L5p0g/gtbbOucag1sBqZh8AayQdG5NGAYsJNa9VNfvjgeficjFwWWwdMBLYGh8ZvAScK6lrrLQ6F3gprtsmaWRsDXBZyr6cc67Rybbn1bXAE7E5y0rg24SgPFXSBGAV8PWYdxqhEXYJUB7zYmabJN1BaB4DcHtVRRZwFfAI0I5QaeUVV865RiurwGpmC4ARaVaNSpPXgKsz7GcyMDlN+lzg+GzK4pxzhc57XjnnXMI8sLrGY9fufJfAuax4YHWNx4ZNtedxrgB4YHWNx8Yt+S6Bc1nxwOoah917YOtH+S6Fc1nxwOoahzK/W3WNhwdW1zhs3AxtW+e7FM5lxQOrK3x7KmDzdujRtfa8zhUAD6yu8JVtATPo2S3fJXEuKx5YXeFbvwnatoGO7fNdEuey4oHVFbbde2DzNujZFTLM2CNptKSlcWqfGzPk+bqkxZIWSfp9Tsvsmj2f/toVtqq2q4enfwwgqQiYBJxDGOJyjqTiOGtAVZ7BwE3A6Wa2WdLhOS61a+b8jtUVtvWboF1b6NAuU45TgRIzWxmnDZpCmB4o1RXAJDPbDGBm63EuhzywusK1azds3Q6HZ34MQOYpf1INAYZI+pukWZLSTlbpM1S4pHhgdYXrw7Lw2qt7fffUkjD775nAJcBvJXWpmclnqHBJ8cDqCpMZfLAROh8WHgVklmnKn1SlQLGZ7TGzd4FlhEDrXE54YHWFadtHsGMXHNGjtpxzgMGSBsYZLsYRpgdK9UfiVO2SehAeDaxMtsDO7eeB1RWmdRuhqEVoZnUQZlYBXEOYU20JMNXMFkm6XdIFMdtLQJmkxcAM4AdmVpbD0rtmzptbucKza3doDdC7BxQV1ZrdzKYR5lpLTbslZdmA78Uf53LO71hd4Vm7Pjxj7dsr3yVx7pB4YHWFZdfuEFh7dq2t0sq5guWB1RUOMyhZE14H9s13aZw7ZFkHVklFkuZL+p/4fqCk2bF/9h9ijSyS2sT3JXH9gJR93BTTl0o6LyW91r7erhlYtS6MuzrgSGjXJt+lce6Q1eWO9TpCrWuVu4D7zGwQsBmYENMnAJtj+n0xH5KGEZrCHAeMBn4Vg3VVX+8xwDDgkpjXNScflsGq90NngH5H5Ls0ztVLVoFVUl/gC8BD8b2As4CnY5ZHgQvj8tj4nrh+VMw/FphiZrtiI+0SQj/vbPp6u6Zsy3ZY+h507ghD+h+s+6pzjUK2d6w/B/4V2Bvfdwe2xDaEUL1/9r6+23H91pg/U5/ubPp6u6aqfCcsKgnjrR53DLTwx/6u8av1Kpb0RWC9mc1rgPLUVhYfJKMp2VMBby8Pd6gnDIJW3qzaNQ3Z3B6cDlwg6T3C1/SzgPuBLpKq/hJS+2fv67sd13cGysjcpzubvt6AD5LR5Ly7FnbuhuMGedMq16TUGljN7CYz62tmAwiVT6+a2TcJXQO/GrONB56Ly8XxPXH9q7HnSzEwLrYaGEgYBOMNsuvr7ZqazdvCICu9uoeBVpxrQurz3esGYIqkfwfmAw/H9IeBxyWVAJsIgZLYf3sqsBioAK42s0oASVV9vYuAyWa2qB7lcoVu63ZYvCI0qRroj9Nd01OnwGpmM4GZcXkloUa/Zp6dwNcybP9j4Mdp0g/o6+2aEDMoWR36/Q/sA8tWhUqq4wZB61b5Lp1zifPaApd76zfB+7Gyccv20BJg6EBo789VXdPkbVtcblVWwoo10CbemW7/OATUDJMDOtcUeGB1ubV2fWhWNeyY/WnDh3onANekeWB1yausDKNULV8VmlR16wydUmr+Wyb7BCrbsSYkXSTJJI1ItADO1eDPWF2yNm6GxStDhRWEaasHHxWWP/3J/ekJSRlr4hxCr705korNbHGNfB0J413MTrQAzqXhd6wuOWbhDrUqePbpBSOOC91VIbQAaNM66aNmO9bEHYQBgXYmXQDnavLA6pJRUQmbtoYa/8FHheeoA49siCPXOtaEpJOBfmb2QkMUyDl/FODqb8dO+PuSEFxbFoXeVFnMVdUQJLUA7gUuzyLvRGAiwFFHHZXbgrkmze9YXf2t2xiCavcucPzghg6qtY010RE4HpgZx7sYCRSnq8DysShcUvyO1dXfxs3QtRMcPygfR9831gQhoI4DvlG10sy2Aj2q3kuaCfyLmc1t4HK6ZsTvWF397N4DO3aFwJoHcczfqrEmlgBT47gUt0u6IC+Fcs2e37G6+tm6PbzmcYSqdGNNmNktGfKe2RBlcs2b37G6+tmyPQyoclj7fJfEuYLhgdXVz9aPwt2qT6ni3D7+1+AO3Z4K+HiHD1TtXA0eWN2h21L1fLVjfsvhXIHxwOoO3QcbQjfVTh3yXRLnCooHVndoduyETdugd09/vupcDf4X4Q7N+xvCmKq9e9Se17lmxgOrq7vKyjDDao8uuRityrlGzwOrq7v3N4SxAY48PN8lca4geWB1dfNReRhztXtnb2blXAa1BlZJ/STNkLRY0iJJ18X0bpKmS1oeX7vGdEl6IE6TsTCOhVm1r/Ex/3JJ41PST5H0VtzmAcknRCpIe/fC0vfC0IDHDvR5q5zLIJs71grg+2Y2jDDk2tWShgE3Aq+Y2WDglfgeYAwwOP5MBH4NIRADtwKnEUZ9v7UqGMc8V6RsN7r+p+YSV/phuGMd0h9a+TATzmVSa2A1s3Vm9ve4vJ0wglAfwvQXj8ZsjwIXxuWxwGMWzAK6SOoNnAdMN7NNZrYZmA6Mjus6mdksMzPgsZR9uUJhFp6tdu0EPbrWnt+5ZqxOz1glDQCGEyZk62Vm6+KqD4BecTnTVBkHSy9Nk+4KyZbtYebVI7x5lXO1yTqwSjoMeAa43sy2pa6Ld5rJTr+ZvgwTJc2VNHfDhg25PpxLtXZ9eLbao0u+S+JcwcsqsEpqRQiqT5jZszH5w/g1nvi6PqZnmirjYOl906QfwKfOyIOPd4SgWrYF+vbyXlbOZSGbVgECHgaWmNm9KauKgaqa/fHAcynpl8XWASOBrfGRwUvAuZK6xkqrc4GX4rptkkbGY12Wsi+Xb/MWQ8lqaNMK+h6R79KkJWm0pKWxVcmNadZ/L7ZqWSjpFUn981FO13xkU7V7OvAt4C1JC2LazcCdwFRJE4BVwNfjumnA+UAJUA58G8DMNkm6gzBHEcDtZrYpLl8FPAK0A16MPy7fyneGSiuIkwQW3t2qpCJgEnAO4fn8HEnFZrY4Jdt8YISZlUu6EvgZcHHDl9Y1F7UGVjP7K5CpweKoNPkNuDrDviYDk9OkzyXMpOkKSdmW8HraCdC2TX7LktmpQImZrQSQNIXQMmVfYDWzGSn5ZwGXNmgJXbNTeLcgrjCYwYbN0KFdIQdVyNzaJJMJ+Dcil2PeytulV7YVtn8Mg5vO40hJlwIjgDMyrJ9I6NTCUUcd1YAlc02N37G69DZvDS0Ajuie75LUJlNrk2oknQ38ELjAzHal25G3OnFJ8cDq0tuyvbFMEjgHGCxpoKTWwDhCy5R9JA0HfkMIquvT7MO5RBX8X43Lg917QouALoU/l5WZVQDXEJrzLQGmmtkiSbdLuiBmuxs4DHhK0gJJxRl251wi/BmrO9DWj8JrIwisAGY2jdDMLzXtlpTlsxu8UK5Z8ztWd6At28MjgMPa57skzjVKHljdgbZsayzPV50rSP6X46or3xl+unbKd0mca7Q8sLrq1sdexod3y285nGvEPLC6/fbuhfVlodLKZ1917pB5YHX7lX4IO3b57KvO1ZMHVhfs3RumXunSEXr61CvO1YcHVhesXR+mXunbq/a8zrmD8sDqoKICVq8LLQG6+9QrztWXB1YHaz6EikoY2Lf2vM65Wnlgbe4+Kg+VVj27QUfvaeVcEnysgKbKLDT0LyoK81UpzSQQm7fBohVh9tWjfcZx55LigbUpMoNFJWGwagiB9cRjoX3b/es/2AjLV0O7NnDCEGjr7VadS4oH1sbOLDwfbVkURvwv3wkfloWBVI48HNq1hhWloUfVgCOhshLeeQ82bobOHeH4Y6ClXwbOJcn/ohqLPRUhcH5Ytj+QVlaGBv3lO8Nd6a49IW/rVnB039B0SoK1G6B8Rwim770PH++ovt45lygPrPlW9Sy0VUtoofB+60ewuyLMjbtjVwig6zeFgAph1KlWLUNwbd0KOh0WKqF6dQ+VUO3bVh+Zqn3bMDHghs3Qri0cP8ibVTmXQwUTWCWNBu4HioCHzOzOPBepbszCz9694c5R7H+t3BuCYmVlSNu9O7xuLw/BdPeezPutuqPs0Bb6HxnGSG3VMlRKZatb53C32+dwOKp3k7tLre3akdQGeAw4BSgDLjaz9xq6nK75KIjAKqkImAScQ5i+eI6kYjNbnHGj8p3w9yXxjYGRZpkQ7GLyvhWW8DaWunEWqmrqu3UKX/G7dgrTTFfuhaIWoUKpY4ew/9Yt9wfCQw2IfQ4PP01QltfOBGCzmQ2SNA64C7i44UvrmouCCKzAqUCJma0EkDQFGAtkDqwifBWuWkaZl5WyUdVyapDKdptqy1QPeFU/LRRGhjILd5ZSCJYt41f3lkU+gHSysrl2xgK3xeWngV9Kklld/yM6l51CCax9gDUp70uB0w66Rbu2cOKQXJbJNQ7ZXDv78phZhaStQHdgY4OU0DU7hRJYsyJpIjAxvt0l6e18lidhPWh6f+i5OKf+Ce9vnwK6vvJ9LTT34x9b3x0USmBdC/RLed83plVjZg8CDwJImmtmIxqmeLnX1M4HGuycsrl2qvKUSmoJdCZUYlVTKNdXvq8FP77m1ncfhfKwbw4wWNJASa2BcYDP/e6ykc21UwyMj8tfBV7156sulwrijjU+97oGeInQZGaymS3Kc7FcI5Dp2pF0OzDXzIqBh4HHJZUAmwjB17mcKYjACmBm04BpddjkwVyVJU+a2vlAA51TumvHzG5JWd4JfK2Ou83n55Hva8GPX0/yb0TOOZesQnnG6pxzTUbBBFZJkyWtT23iIukkSbMkLZA0V9KpMV2SHpBUImmhpJNTthkvaXn8GZ/uWA0lwzl9UtLrkt6S9LykTinrborntFTSeSnpo2NaiaQbG/o8UsrRT9IMSYslLZJ0XUzvJml6/J1Pl9Q1phfc51Tb71JSG0l/iOtnSxqQsi7t55Pw8b8Xf78LJb0iqX/Kusr4t7BA0iFV7mZx/MslbUg5zj+mrKvXZ5bFse9LOe4ySVtS1iVx7gf8PdZYn9z1amYF8QN8HjgZeDsl7c/AmLh8PjAzZflFQj+okcDsmN4NWBlfu8blrgV2TnOAM+Lyd4A74vIw4E2gDTAQWEGojCmKy0cDrWOeYQc55s2E/vK5OJ/ewMlxuSOwLJb7Z8CNMf1G4K5C/Jyy+V0CVwH/GZfHAX842OeTg+P/A9A+Ll9Zdfz4/qMGOP/LgV+m2bZen9khXMfXEioiEzn3uI8D/h5rrE/sei2YO1Yze41QY1stGai6o+sMvB+XxwKPWTAL6CKpN3AeMN3MNpnZZmA6MDr3pU8vwzkNAV6Ly9OBi+LyWGCKme0ys3eBEkJ3zX1dNs1sN1DVZTPTMX9iZv+YaX19mNk6M/t7XN4OLCH0ahoLPBqzPQpcmHJOhfQ5ZfO7TD2Xp4FRkkTmzyfR45vZDDMrj29nEdrlJqVO11IN9f3M6nrsS4An67D/WmX4e0yV2PVaMIE1g+uBuyWtAe4Bborp6box9jlIeiFZxP4L6mL2N25vVOcUvyIPB2YDvcxsXVz1AVA1h3ahnVM2x63W/RWo6v6aRJnruo8JhDuoKm3jI7FZki7MtFECx78ofhV+WlJt12fSxyY+/hgIvJqSXN9zr08Z63zuhR5YrwS+a2b9gO8S2iM2dv2Bn0kqB+4EDpM0qGqlpEck/Xt825Hw1fpbkjZJ+l+qhomRbpC0VtL2+NxqVEy/TdLvUvb3lKQPJG2V9Jqk4+p7ApIOA54BrjezbanrLHx38qYm9STpUmAEcHdKcn8LPZK+Afxc0jE5OPTzwAAzO5FwZ/ZoLflzYRzwtJlVpqQ1xLknptAD63jg2bj8FPu/emXqxphV19g82wOUE/ojD09JT1f2Mwn/Kf9CuAu8mfCfchdwDfApM+tI+KryXobjvQgMBg4H/g48UZ/CS2pFCKpPmFnVZ/Nh/MpEfF1/kHPK5+dUl+6vqHr31yTKnNU+JJ0N/BC4wMx2VaWb2dr4uhKYSfXrJ5Hjm1lZyjEfIoxhm3XZ63PsFOOo8RgggXPPRnLXa30fCCf5AwygekXPEuDMuDwKmBeXv0D1h8xvpDxkfpfwgLlrXO5WYOe0mlBp1YIw+LIBg4DjCA/0HwN+QXhAfgehO+Zqwlejqof+ownB62ygVY3j3Qb8LkNZusTjdT7Ec1Es389rpN9N9cqrnxXi50ToELOyxu/yuBp5rqZ65dXUuFz1+VRVXq2k7pVX2Rx/OKGSZ3CN9K5Am7jcA1jOQSp/6nH83inLXwZmJfGZZXPsmG8o4SZBSZ57pr/HGusSu15zeiHX8YSfBNYR7uhKCc+XPgvMix/CbOCUmFeEwY1XAG8BI1L28x1CxUIJ8O0CPKeyuLyM8CjAgEEx/w+BbTHPGMKjgP+I+9gd038Y834D+CuwmVARcGRMv40YWAk1sXfG39M2YEs83jGHeD6fjdsvBBbEn/MJzyBfiRf8y1UXXSF+TrG8y2KZqn6XtxPuDgHaEr4dlQBvAEenbPvDuN1SYmuVHBz/ZeDDlN9vcUz/TPwdvhlfJ+To+D8l1AO8CcwAhib1mdV27JTr984a2yV17un+Hv8J+Kekr1fvedXAJL0H/KOZvRzffwx82swWxvd/IvRx/1GN7Y4nPMy/xMxeSUnvBPwGqDCzb0m6jRCoL5X0LcLjg/MJdwGdCYF4sJmV5PREnWvGCv0Za3OwAPiGpCKFuZvOqFoh6YuSBsXmPluBSmCvpGMlnaUwl9NOYAewN82+OxKex5YB7YGf5PhcnHN4YC0E1wFfInxN/ybwx5R1gwlfDT8CXgd+ZWYzCM/57iQMBvwBoWLqJg70GLCK8KB9MaFdpHMux/xRgHPOJczvWJ1zLmFZBVZJ31UYdONtSU9KaqswYvvsOGDBHxRGbz+kQSxqG5zBuUzqM7CGc7lSa2CV1Af4Z0LTg+MJTXiq5ma/z8wGEWqaJ8RN9s3hDtwX8yFpWNzuOEI7zF/FCpuqeeHHEAa6uCTmdS4bj3DwfttjCM+qBxMmCvx1A5TJNXPZPgpoCbSLPVHaE9qCnUUYpAIOHHijLoNY1GdgCNfM2aEPrOFcztQaWC10JbuH0PtnHaHZzzxgi4VBKqD6oAR1HcSiIAcZcU2GX1+uwdU655XCoMVjCV3RthB6peRlKD6lzPveoUOHU4YOHZqPYrg8mjdv3kYz65mLffv15SCZayybyQTPBt41sw0Akp4FTid8pWoZ70pTByXINIf7wQYyyGqAA0uZ933EiBE2d269p/92jYykVXXcJOsBNPz6cnBI19gBsnnGuhoYKal9fFY6itDYfAZhjnYIo1A9F5czzeFeDIyLrQYGEioT3iC7eeGdO1TFwGWxdcBIYKvtHzvWuZyo9Y7VzGZLepow5FwFMJ/wX/0FYEocO3Q++8dKTTuHu4W53qcSgnIFcLXF8RaVZl745E7RNWWSniQMr9hDUilwK9AKwMz+kzAt9vmEytJy4Nv5KalrThptzyv/qtY8SZpnYcDjnPLrq/lK4hrznlfOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6zOOZcwD6yu0ZM0WtJSSSWSbkyz/ihJMyTNl7RQ0vn5KKdrPjywukZNUhEwCRgDDAMukTSsRrYfAVPNbDhhqqBfNWwpXXPjgdU1dqcCJWa20sx2A1MI07WnMqBTXO4MvN+A5XPNUDbTXztXyPoAa1LelwKn1chzG/BnSdcCHQhTujuXM37H6pqDS4BHzKwvYcbWxyUdcO1LmihprqS5GzZsaPBCuqbDA6tr7NYC/VLe941pqSYAUwHM7HWgLdCj5o7M7EEzG2FmI3r27Jmj4rrmIKvAKqmLpKclvSNpiaRPS+omabqk5fG1a8wrSQ/EGtqFkk5O2c/4mH+5pPEp6adIeitu84AkJX+qromaAwyWNFBSa0LlVHGNPKuBUQCSPkEIrH5L6nIm2zvW+4E/mdlQ4JPAEuBG4BUzGwy8Et9DqJ0dHH8mAr8GkNQNuJXw/OtU4NaqYBzzXJGy3ej6nZZrLsysArgGeIlwXU41s0WSbpd0Qcz2feAKSW8CTwKXm5nlp8SuOai18kpSZ+DzwOUAseZ1t6SxwJkx26PATOAGQo3sY/HCnRXvdnvHvNPNbFPc73RsA1IkAAAWwElEQVRgtKSZQCczmxXTHwMuBF5M5Axdk2dm04BpNdJuSVleDJze0OVyzVc2d6wDCV+b/is2sH5IUgegl5mti3k+AHrF5XS1tH1qSS9Nk+6cc41SNoG1JXAy8OvYwPpj9n/tByDeneb8q5XX2jrnGoNsAmspUGpms+P7pwmB9sP4FZ/4uj6uz1RLe7D0vmnSD+C1ts65xqDWwGpmHwBrJB0bk0YBiwk1r1U1++OB5+JyMXBZbB0wEtgaHxm8BJwrqWustDoXeCmu2yZpZGwNcFnKvpxzrtHJtufVtcATsTnLSuDbhKA8VdIEYBXw9Zh3GqERdglQHvNiZpsk3UFoHgNwe1VFFnAV8AjQjlBp5RVXzrlGK6vAamYLgBFpVo1Kk9eAqzPsZzIwOU36XOD4bMrinHOFznteOedcwjywOudcwjywOudcwjywOudcwjywOudcwjywOudcwjywOudcwjywOudcwjywOudcwjywOudcwjywukZP0mhJS+PUPjdmyPN1SYslLZL0+4Yuo2tefPpr16hJKgImAecQhricI6k4zhpQlWcwcBNwupltlnR4fkrrmgu/Y3WN3alAiZmtjNMGTSFMD5TqCmCSmW0GMLP1OJdDHlhdY5dpyp9UQ4Ahkv4maZYkn6zS5ZQ/CnDNQUvC7L9nEmaoeE3SCWa2JTWTpImEmYU56qijGrqMrgnxO1bX2GWa8idVKVBsZnvM7F1gGSHQVuNT/7ikeGB1jd0cYLCkgXGGi3GE6YFS/ZE4VbukHoRHAysbspCuefHA6ho1M6sAriHMqbYEmGpmiyTdLumCmO0loEzSYmAG8AMzK8tPiV1z4M9YXaNnZtMIc62lpt2SsmzA9+KPcznnd6zOOZcwD6zOOZcwD6zOOZewrAOrpCJJ8yX9T3w/UNLs2D/7D7FGFklt4vuSuH5Ayj5uiulLJZ2Xkl5rX2/nnGss6nLHeh2h1rXKXcB9ZjYI2AxMiOkTgM0x/b6YD0nDCE1hjgNGA7+Kwbqqr/cYYBhwSczrnHONUlaBVVJf4AvAQ/G9gLOAp2OWR4EL4/LY+J64flTMPxaYYma7YiPtEkI/72z6ejvnXKOR7R3rz4F/BfbG992BLbENIVTvn72v73ZcvzXmz9SnO5u+3s4512jUGlglfRFYb2bzGqA8tZVloqS5kuZu2LAh38VxDWxZ2bJ8F8G5rGRzx3o6cIGk9whf088C7ge6SKrqYJDaP3tf3+24vjNQRuY+3dn09Qa8L3dzd+/r9+a7CM5lpdbAamY3mVlfMxtAqHx61cy+Sega+NWYbTzwXFwuju+J61+NPV+KgXGx1cBAwiAYb5BdX2/XzK3/eD2PLHgk38VwLiv16dJ6AzBF0r8D84GHY/rDwOOSSoBNhEBJ7L89FVgMVABXm1klgKSqvt5FwGQzW1SPcrkmaNIbk9hVuSvfxXAuK3UKrGY2E5gZl1cSavRr5tkJfC3D9j8Gfpwm/YC+3s5VKd9TzqQ5k7jg2Aso9i8zrhHwnleu4D2y4BHKdpTxg8/8IN9FcS4rHlhdQavYW8G9r9/LaX1O4/R+p+e7OM5lxQOrK2i/f+v3rNi8ghtOv4HQz8S5wueB1RWsir0V3PHaHZx0xElcOPTC2jdwrkB4YHUF64mFT1CyqYRbz7j1oHer2Q7iI+kiSSZpRE4K7FzkgdUVpKq71eFHDGfssZmHjsh2EB9JHQkDCc3OUZGd28cDqytIj7/5OCs2r+C2M2+r7dlqtoP43EEYaW1n8qV1rjoPrK7glO8p55aZtzDiyBF8aciXaste6yA+kk4G+pnZCwfbkY9F4ZLigdUVnHv+3z2Ubivl3nPvrXdLAEktgHuB79eW18eicEnxwOoKSum2Uu762118bdjX+Fz/z2WzSW2D+HQEjgdmxoGERgLFXoHlcskDqysoN71yE5V7K/nZOT/LdpODDuJjZlvNrIeZDYgDCc0CLjCzuUmX3bkqHlhdwZhdOpvfLfwd3//09xnQZUBW28TB1KsG8VkCTI0D/twu6YLclda5zOozupVziTEzrn/peo447Ahu/Gzd5pNMN4iPmd2SIe+Zh1xI57LkgdUVhClvT2FW6SwmXzCZjm065rs4ztWLPwpweVe+p5wbXr6B4UcMZ/xJ42vfwLkC53esLu9+8r8/Yc22NfzuK7+jhfx/vWv8/Cp2De7dze/yv6v+F4A31r7BXX+7i8s+eRmf7//5PJfMuWT4HatrUFt2buFTv/0UZTvK+OWYX/KjGT+i92G9ufdcnyjQNR0eWF2DenTBo5TtKAPgmhevoXVRa16f8Drd23fPc8mcS44/CnANYu77c9G/iZtfvZkRR45g5viZjOw7kkcvfJShPYbmu3jOJcrvWF2D+O283wKhBcBXP/FVzhhwBq9PeD3PpXIuN/yO1TWI9eXr9y37bACuqas1sErqJ2mGpMWSFkm6LqZ3kzRd0vL42jWmS9IDcTT3hXHItqp9jY/5l0san5J+iqS34jYPyCc3alLMjNml+8eXPrbHsXksjXO5l80dawXwfTMbRhgZ6Oo4QvuNwCtmNhh4Jb6HMJL74PgzEfg1hEAM3AqcRhic+NaqYBzzXJGy3ej6n5rLtT+v+DM3TL+BS565hI3lGzPmK91WyrqP1nHdadex7JplDVhC5/Kj1mesZrYOWBeXt0taQhhIeCxwZsz2KDATuCGmP2ZmBsyS1EVS75h3upltApA0HRgtaSbQycxmxfTHgAuBF5M5RZekj3Z/xAvLXuCZJc/w1OKn9qX3aNeDX5z/i7TbzCqdBcClJ17K4O6DG6SczuVTnZ6xShoADCfMG9QrBl2AD4BecTnTiO4HSy9Nk+4KzJqtaxj0wCDGPTOOP77zR7467KvM/z/z+cYJ3+DJt5+kYm9F2u1mlc6ibcu2nNjrxAYusXP5kXVglXQY8AxwvZltS10X704t4bKlK4NPnZEn/7Psfxg6aSibd27myYueZMuNW3jqa09x0hEn8ZWhX6FsRxlvrH0j7bavl77OKb1PoXVR6wYutXP5kVVgldSKEFSfMLNnY/KH8Ss+8bWq2jfTiO4HS++bJv0APnVGfpRuK+Wbz36Twd0GM2P8DMYdP472rdrvW39CrxMAWLFpxQHbbtu1jTnvz+FzR2U1G4BzTUI2rQIEPAwsMbPUfofFQFXN/njguZT0y2LrgJHA1vjI4CXgXEldY6XVucBLcd02SSPjsS5L2ZfLMzPjiuevoGJvBc9e/Cyf6feZA/L06xT+X67auuqAdTPenUHF3grOG3RezsooabSkpbFVyQGDuUr6XmzVslDSK5L656wwzpFdB4HTgW8Bb0laENNuBu4EpkqaAKwCvh7XTQPOB0qAcuDbAGa2SdIdhKk0AG6vqsgCrgIeAdoRKq284qpATJ4/mT+V/IlfjPkFR3c9Om2edq3acXiHw1m15cDA+ucVf6ZDqw5pA3ISJBUBk4BzCM/n50gqNrPFKdnmAyPMrFzSlcDPgItzUiDnyK5VwF+BTO1KR6XJb8DVGfY1GZicJn0uYcI3V0DMjDv/dicj+47kqk9dddC8/Tv3T3vH+tKKl/iHgf+Qy+erpwIlZrYSQNIUQsuUfYHVzGak5J8FXJqrwjgH3vPKHcTb69+mZFMJl55waa3jpA7rOYz5H8wn/F8NVmxawYrNKzjvmNw9BiBza5NMJuDfiFyOeWB1ae2p3MOl/30pXdp24cuf+HKt+T931OfYWL6Rdza+sy/t+WXPA+Q6sGZN0qXACODuDOu91YlLhAdWl9aUt6ew8MOF/NfY/+LIjkfWmr/qGerstaHr6l7by6Q5kzitz2m57hSQqbVJNZLOBn5ImPp6V7odeasTlxQPrC6tKYum0KdjH8YeOzar/EO6D6Fty7a89eFbALyw7AVKNpXw3ZHfzWUxIVSGDpY0UFJrYByhZco+koYDvyEE1fVp9uFcojywugM8/ubjTFs+jUtPvJRsx8MpalHEcT2PY8GHoeHIs+88S7d23fjKJ76Sy6JiZhXANYTmfEuAqWa2SNLtki6I2e4GDgOekrRAUnGG3TmXCB+P1VVTubeSO167g+FHDOfHZ/24TtueffTZ3PW3u/i3mf/GC8te4KyBZ9GqqFWOSrqfmU0jNPNLTbslZfnsnBfCuRR+x+qqeXrx0yzftJybP3czRS2K6rTtjz7/I/p37s9tf7mNDeUbGDNoTI5K6Vxh88Dq9tlduZubX72Z4w8/ni8Prb0lQE2HtT6Mly97mZ7te9K+VXvGHT8uB6V0rvD5owC3z6Q3JrFy80r+9M0/1flutcqgboNY9/11lO8przaegHPNid+xOgDmr5vPbX+5jfOOOa/e/fqLWhTRsU3HhErmXOPjgdWxu3I3lzxzCZ3adOJXX/hVvovjXKPnjwKauVVbVnHzqzeztGwp074xLeNAK8657HlgbcZ++cYvufbFawH48tAvM2aw1+I7lwQPrM3Uyytf5toXr6Vbu26MHjSae8+9t/aNnHNZ8cDazBQvLeaRBY/w/LLnGdpjKH+f+HfatWqX72I516R4YG1G/vLeXxg7JfT9v3LElfzfz/9fD6rO5YAH1iZoVuksnlr0FF3aduH6kdfTsU1H9lTu4coXrqRvp768dvlrDOw6MN/FdK7J8sDaxBQvLeaiqRftm4r6p3/9KV869kuUlZexZOMSiscVe1B1Lsc8sDYh67av49JnL+WkI07ivy/+b9ZsXcPvFv6OpxY/xYbyDXzlE1/hS8d+Kd/FdK7J88BaIPbaXuavm8/qrasZM3gMbYra8NfVf2Xrrq3MLp3NkR2PZOIpEw/oajrv/Xk8OO9BZrw3g7IdZVTsreCJrzxB30596dupL5/u92nuH3M/s0pnMfyI4Xk6O+eaFw+sebJu+zqeWfIMO/bsYNbaWSxav4ilZUsB+NSRn6J8TzmLNiyqts3yTcu597z9zaKu/9P13D/7ftq1bMc5x5xDyxYt+d7I7zGk+5Bq27Vs0ZLPHvXZ3J+Ucw7wwFpvFXsr2Fmxkx17drCzYifle8r3/Xy85+PwuvtjVm9dzV9W/YUBXQbwiR6f4D9e/w/WbAtz4HVt25UTep3AkO5D+HjPx7y9/m2KVMQVJ1/BoG6D+NaJ3+K2mbdx36z7+NaJ32J47+E8985z3D/7fiaePJG7zrmLLm275Pk34ZyrUjCBVdJo4H6gCHjIzO48WP4dFTtY8MECKvdWUmmV1V4r9lYckFZpMX1vJbsrd7OjYse+gJh2uWLHvmB5sOWqSqJsHdb6MD7a/RHd23Xn1jNuZdzx4xjSfUi1WVAr91bSQi2qjd7/41E/5qH5D3HP6/dw4bEXcuULV3Jy75P5xfm/yOXU0o1CbdeOpDbAY8ApQBlwsZm919DldM1HQQRWSUXAJOAcwvTFcyQVm9niTNssXr+Y4b9J7plhu5btaNuyLe1atTtguWPrjhze4fD96S3bpc1XtdyhVQfat2pPh9bxtVUH9tpejux4JG1atmHmezM56YiTOLzD4WnLkm7Ivh7te3DxcRfz+7d+z+/f+j39O/fnyYue9KCa3bUzAdhsZoMkjQPuAi5u+NK65qIgAitwKlBiZisBJE0BxgIZA+vR3Y7mnq/fQ1GLIopUtO+1ZYuWB6TVXNe6qPW+4Ni2ZVvaFLXJem6nJJx7zLmHtN3jX36cqz51Fau3ruaLQ75IpzadEi5Zo5TNtTMWuC0uPw38UpLMzBqyoK75KJTA2gdYk/K+FDjtYBt0bds1q/num5KiFkVeCXWgbK6dfXnMrELSVqA7sLFBSuianUIJrFmRNBGYGN/ukvR2PsuTsB40vT/0XJxT/4T3t08BXV/5vhaa+/GPre8OCiWwrgX6pbzvG9OqMbMHgQcBJM01sxENU7zca2rnAw12TtlcO1V5SiW1BDoTKrGqKZTrK9/Xgh9fc+u7j0KZQWAOMFjSQEmtgXGAz/3uspHNtVMMjI/LXwVe9eerLpcK4o41Pve6BniJ0GRmspktqmUz5zJeO5JuB+aaWTHwMPC4pBJgEyH4OpczBRFYAcxsGjCtDps8mKuy5ElTOx9ooHNKd+2Y2S0pyzuBr9Vxt/n8PPJ9Lfjx60n+jcg555JVKM9YnXOuySiYwCppsqT1qU1cJJ0kaZakBZLmSjo1pkvSA5JKJC2UdHLKNuMlLY8/49Mdq6FkOKdPSnpd0luSnpfUKWXdTfGclko6LyV9dEwrkXRjQ59HSjn6SZohabGkRZKui+ndJE2Pv/PpkrrG9IL7nGr7XUpqI+kPcf1sSQNS1qX9fBI+/vfi73ehpFck9U9ZVxn/FhZIOqTK3SyOf7mkDSnH+ceUdfX6zLI49n0px10maUvKuiTO/YC/xxrrk7tezawgfoDPAycDb6ek/RkYE5fPB2amLL8ICBgJzI7p3YCV8bVrXO5aYOc0BzgjLn8HuCMuDwPeBNoAA4EVhMqYorh8NNA65hmWp/PpDZwclzsCy2K5fwbcGNNvBO4qxM8pm98lcBXwn3F5HPCHg30+OTj+PwDt4/KVVceP7z9qgPO/HPhlmm3r9ZnV9ToGriVURCZy7nEfB/w91lif2PVaMHesZvYaoca2WjJQdUfXGXg/Lo8FHrNgFtBFUm/gPGC6mW0ys83AdGB07kufXoZzGgK8FpenAxfF5bHAFDPbZWbvAiWE7pr7umya2W6gqstmgzOzdWb297i8HVhC6NU0Fng0ZnsUuDAuF9rnlM3vMvVcngZGSRKZP59Ej29mM8ysPL6dRWiXm5T6XEv1/czqeuxLgCfrsP9aZfh7TJXY9VowgTWD64G7Ja0B7gFuiunpujH2OUh6IVnE/gvqa+xv3N6ozil+RR4OzAZ6mdm6uOoDoFdcLrRzyua41bq/AlXdX5Moc133MYFwB1WlbXwkNkvShZk2SuD4F8Wvwk9Lqu36TPrYxMcfA4FXU5Lre+71KWOdz73QA+uVwHfNrB/wXUJ7xMbuO8BVkuYRvk7vznN56kzSYcAzwPVmti11nYXvTt7UpJ4kXQqMAO5OSe5voUfSN4CfSzomB4d+HhhgZicS7swerSV/LowDnjazypS0hjj3xBR6YB0PPBuXn2L/V69M3Riz6hqbT2b2jpmda2anEL7qrIirGsU5SWpFCKpPmFnVZ/Nh/MpEfF0f0wvtnOrS/RVV7/6aRJmz2oeks4EfAheY2a6qdDNbG19XAjMJ3xgSPb6ZlaUc8yHCGLZZl70+x04xjhqPARI492wkd73W94Fwkj/AAKpX9CwBzozLo4B5cfkLVH/I/EbKQ+Z3CQ+Yu8blbgV2TofH1xaEwZe/E98fR/XKkZWEB/4t4/JA9j/0Py5P56JY5p/XSL+b6pVXPyvEzymb3yVwNdUrr6Ye7PPJwfGHE/7ZDq6R3hVoE5d7AMupYyVmlsfvnbL8ZWBWEp9ZttcxMBR4j9jGPqlzT9lXtb/HGusSu14b/I/zICf8JLAO2EN4hjEB+CwwL34Is4FTYl4RBjdeAbwFjEjZz3cIFQslwLcL8JyuI9SmLwPurHEB/TCe01Jia4iYfn7MvwL4YR7P57OEr/kLgQXx53zCM8hX4gX/ctVFV4ifU7rfJXA74e4QoC3h21EJ8AZwdG2fT8LHfxn4MOX3WxzTPxN/h2/G1wk5Ov5PCfUAbwIzgKFJfWa1HTu+vw24s8Z2SZ17ur/HfwL+Kenr1XteOedcwgr9GatzzjU6Hlidcy5hHlidcy5hHlidcy5hHlidcy5hHlidcy5hHlidcy5hHlidcy5h/x8PJYRu0c/wpQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 4 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\r\n",
    "import pandas as pd\r\n",
    "data=pd.read_csv(\"line_animation.csv\")\r\n",
    "\r\n",
    "year=data['time']\r\n",
    "japan=data['japan']\r\n",
    "russia=data['russia']\r\n",
    "data.describe()\r\n",
    "\r\n",
    "fig2,axes2=plt.subplots(2,2)\r\n",
    "fig2.set_size_inches(5,5)\r\n",
    "ax1=axes2[0,0]\r\n",
    "ax1.set_title(\"japan\")\r\n",
    "ax1.set_xlim([1800,2040])\r\n",
    "ax1.set_ylim([0,80000])\r\n",
    "ax1.plot(year,japan,color='pink')\r\n",
    "ax2=axes2[1,0]\r\n",
    "ax2.set_title(\"russia\")\r\n",
    "ax2.set_xlim([1800,2040])\r\n",
    "ax2.set_ylim([0,80000])\r\n",
    "ax2.plot(year,russia,color='green')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "py35-paddle1.2.0"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
